Dispersive wave runup on non-uniform shores 
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Abstract Historically the finite volume methods have been developed for the numer- 
ical integration of conservation laws. In this study we present some recent results on 
the application of such schemes to dispersive PDEs. Namely, we solve numerically 
a representative of Boussinesq type equations in view of important applications to 
the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto 
a non-uniform beach are presented along with great lines of the employed numerical 
method (see D. Dutykh et al. (2011) ||6l for more details). 
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1 Introduction 

The simulation of water waves in realistic and complex environments is a very chal- 
lenging problem. Most of the applications arise from the areas of coastal and naval 
engineering, but also from natural hazards assessment. These applications may re- 
quire the computation of the wave generation 15] [121, propagation ifTTI . interaction 
with soUd bodies, the computation of long wave runup lfT6l [TSll and even the ex- 
traction of the wave energy ifTSl . Issues like wave breaking, robustness of the nu- 
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merical algorithm in wet-dry processes along with the validity of the mathematical 
models in the near-shore zone are some basic problems in this direction fTT}. Dur- 
ing past several decades the classical Nonlinear Shallow Water Equations (NSWE) 
have been essentially employed to face these problems Q. Mathematically, these 
equations represent a system of conservation laws describing the propagation of in- 
finitely long waves with a hydrostatic pressure assumption. The wave breaking phe- 
nomenon is commonly assimilated to the formation of shock waves (or hydraulic 
jumps) which is a common feature of hyperbolic PDEs. Consequently, the finite 
volume (FV) method has become the method of choice for these problems due to its 
excellent intrinsic conservative and shock-capturing properties l2\Ul- 

In the present article we report on recent results concerning the extension of the 
finite volume method to dispersive wave equations steming essentially from water 
wave modehng lfT4l l4ll6l. 



2 Mathematical model and numerical methods 

Consider a cartesian coordinate system in two space dimensions {x,z) to simplify 
notations. The z-axis is taken vertically upwards and the x-axis is horizontal and 
coincides traditionally with the still water level. The fluid domain is bounded below 
by the bottom z = —h{x) and above by the free surface z = r\{x,t). Below we will 
also need the total water depth H{x,t) := h{x) + I7(x,f). The flow is supposed to 
be incompressible and the fluid is inviscid. An additional assumption of the flow 
irrotationality is made as well. 

In the pioneering work of D.H. Peregrine (1967) lfT4l the following system of 
Boussinesq type equations has been derived: 

^, + {{h + ^)u)^ = Q, (1) 

h 

u, + M«v + gr\x - T {hu)xxt + -rUxxt = 0, (2) 
2 O 

where u{x,t) is the depth averaged fluid velocity, g is the gravity acceleration and 
underscripts (m^, r\t) denote partial derivatives. 

In our recent study |l6l we proposed an improved version of this system which 
contains higher order nonlinear terms which should be neglected from asymptotic 
point of view and can be written in conservative variables {H,Q) — {H,Hu) as: 

H, + Q, = 0, (3) 

((1 + - ^HH,,)Qr - \^H^Q,^ - \^HH,Q^) + + f , " 

Obviously the linear characteristics of both systems (O and (O, dUi coincide 
since they differ only by nonlinear terms. 
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However, this modification has several important impUcations onto structural 
properties of the obtained system. First of all, the magnitude of the dispersive terms 
tends to zero when we approach the shoreline H ^0. This property corresponds to 
our physical representation of the wave shoaling and runup process. On the other 
hand, the resulting system becomes invariant under vertical translations (subgroup 
Gs in Theorem 4.2, T. Benjamin & P. Olver (1982) 0): 

z^z + d, 77-s— T]— c/, h^h + d, iii-u, (5) 

where d is some constant. This property is straightforward to check since we use 
only the total water depth variable H = h + r] which remains invariant under trans- 
formation (|5]). 

Remark 1 In this paper we will consider the initial-boundary value problem posed 
in a bounded domain I = [b\^b2] with reflective boundary conditions. In this case 
one needs to impose boundary conditions only in one of the two dependent variables, 
cf. SiEil- In the case of reflective boundary conditions it is sufficient to take u(bi , t) = 
u{b2,t)=Qi. 



2.1 Finite volume discretization 

Let .'7 = {xi}, i E Z denotes a partition of M into cells C,- = 1 ,x.^i ) where 
Xi = (x. , i +x. i)/2 denotes the midpoint of d . Let Ax, =x. , i —x. 1 be the length 
of the cell C,, Ax.^ 1 = .y,+i — x,. (Here, we consider only uniform grids with Axj = 
Ax.^^ = Ax.) 

The governing equations (O, (01 can be recast in the following vector form: 



where 



[D(yt)] + [F(y)], = S(y), 



- 1,(1 + - IHHM \H^Q.. - IhHm J ' ^"'^ 

We denote by H, and U, the corresponding cell averages. To discretize the dis- 
persive terms in ^ we consider the following approximations; 
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Ax Jx. 



i o 



Qdx: 



1 

3 V 2Ax 



1^ Hi+i-2Hi + Hi^i \ ^ 
-6^' A? P'-' 



- r 

AxJx. 



^ ^HHxQx dx 



1^^ H,+i-H,^ 
3 ' 2Ax 2Ax 



Qi+i - Qi- 



1 

Ax 



1 1 T 



1 2 a+i-2a+ a 

3 ' 4x2 



We note that we approximate the reflective boundary conditions by taking the 
cell averages of u on the first and the last cell to be mq = "iV+i = 0. We do not 
impose explicitly boundary conditions on H. The reconstructed values on the first 
and the last cell are computed using neighboring ghost cells and taking odd and even 
extrapolation for u and H respectively. These specific boundary conditions appeared 
to reflect incident waves on the boundaries while conserving the mass. 

This discretization leads to a linear system with tridiagonal matrix denoted by L 
that can be inverted efficiently by a variation of Gauss elimination for tridiagonal 
systems with computational complexity 0{n), «-being the dimension of the system. 
We note that on the dry cells the matrix becomes diagonal since //, is zero on dry 
cells. For the time integration the explicit third-order TVD-RK method is used. In 
the numerical experiments we observed that the fully discrete scheme is stable and 
preserves the positivity of H during the runup under a mild restriction on the time 
step At. 

Therefore, the semidiscrete problem of (|6]l 
in the form; 



(|7]i is written as a system of ODEs 
1 



Ax' 

where L,- is the th row of matrix L and i 



l) 



-Si, 



3' Ax 

can be chosen as one of the numer- 



ical flux functions |6 1 (in computations presented below we choose the FVCF flux 
191). In the sequel we will use the KT and the CF numerical fluxes. In this case the 
Jacobian of F is given by the matrix 



1 

gH-{Q/Hf 2Q/H^ 



and the eigenvalues are X\^2 — Q/H±^/gH. Therefore, the characteristic numerical 
flux 191 takes the form 



F(V^ 



■-u(M)- 



-F(V[ 



where ii — {^i^iizY are the Roe average values. 
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2 



and 



U(M) 




) 



c 



'g^ir, ii = sign(A,). 



For more details on the discretization and reconstruction procedures, (that are 
based on the hydrostatic reconstraction, IT]), we refer to our complete work on this 
subject 161. 



3 Numerical results 

In the present section we show a numerical simulation of a solitary wave runup onto 
a non-uniform sloping beach. More precisely, we add a small pond along the slope. 
As our results indicate, this small complication is already sufficient to develop some 
instabilities which remain controlled in our simulations. 

As an initial condition we used an approximate solitary wave solution of the 
following form: 



where is the amplitude relative to the constant water depth taken to be unity in 
our study. The solitary wave speed along with the wavelength X are given here: 



The solitary wave is centered initially at xq = 10.62 and has amplitude A, = 0.08. 
The constant slope j5 is equal to 2.88°. The sketch of the computational domain can 
be found in [61. 

In numerical simulations presented below we used a uniform space discretization 
with Ax = 0.025 and very fine time step At ^ Ax / lOQ to guarantee the accuracy and 
stability during the whole simulation. 

Snapshots of numerical results are presented on Figures [T]-|6] We present simul- 
taneously three different computational results: 

• Modified Peregrine system solved with UN02 reconstruction ifTol 

• The same system with classical MUSCL TVD2 scheme lfT3l 

• NonUnear Shallow Water Equations (NSWE) with UN02 scheme ifTOl 



r]o{x) =AsSech^{X{x~xo)), uo{x) 



1 + rioix) ' 
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(a) f = 1 s (b) ; = 3 s 

Fig. 1 Solitary wave aproaching a sloping beach with a pond. 
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(a) r-3.5s (b) f-4s 

Fig. 2 Beginning of the pond inundation. 
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X X 



{a)f = 5s (b)f = 5.5s 

Fig. 3 A part of the wave mass is trapped in the pond volume. 

Surprisingly good agreement was obtained among all three numerical models. 
Presumably, the complex runup process under consideration is governed essen- 
tially by nonlinearity. However, on Figures [Ttb) and |2ja) the amplitude predicted 
by NSWE is slightly overestimated. 

On Figures|3lb) -SJb) some oscillations (due to the small-dispersion effect char- 
acterizing dispersive wave breaking procedures) can be observed. However, their 
amphtude remains small for all times and does not produce any blow up phenom- 
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— Boussinesq (UNO) 

- Boussinesq (TVD) 
-Shallow water (UNO) 



(a) f = 6 s 
Fig. 4 Wave oscillations in the pond. 




(b) f = 6.5 s 




-Boussinesq (UNO) 
- Boussinesq (TVD) 
-Shallow water (UNO)[ 




-2 2 A 6 B 10 -2 

(a) t = l s 
Fig. 5 Stabilization of wave oscillations. 


2 

(b) f = 


0.3 




Boussinesq (UNO) 

- - Boussinesq (TVD) 
Shallow water (UNO) 


0.2 
0.1 







-0.1 
-0.2 




-2 2 4 
X 

Fig. 6 The whole system is tending to the rest position (t 


6 8 1 

= 10 s). 



ena. Later these oscillations decay tending gradually to the "lake at the rest" state 
(see Figures|5]|6]l. 

In the specific experiment a friction term could be beneficial to reduce the ampli- 
tude of oscillations (or damp them out completely). However, we prefer to present 
the computational results of our model without adding any ad-hoc term to show its 
original performance. 
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In this study we presented an improved version of the Peregrine system which is par- 
ticularly suited for the simulation of dispersive waves runup. This system allows for 
the description of higher amplitude waves due to improved nonlinear characteristics. 
Better numerical stability properties have been obtained since most of the disper- 
sive terms tend to zero when we approach the shoreline. Consequently, our model 
naturally degenerates to classical Nonlinear Shallow Water Equations (NSWE) for 
which the runup simulation technology is completely mastered nowadays. However 
we underline that there is no artificial parameter to turn off dispersive terms. Their 
importance is naturally governed by the underlying physical process. 
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